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Abstract 



We derive a closed-form (infinite series) representation for the dis- 
tribution of the ruin time for the Sparre Andersen model with expo- 
nentially distributed claims. This extends a recent result of Dickson et 
al. [7] for such processes with Erlang inter-claim times. We illustrate 
our result in the cases of gamma and mixed exponential inter-claim 
time distributions. 
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1 Introduction 

In the Sparre Anderson model, the (continuous-time) surplus process {U(t)} t >o 
has the form 



where u > is the initial surplus, c > is the premium rate, and {N(t)} t >o 
is a delayed renewal process generated by a sequence of inter-claim times 



and {Xj}j>i is the sequence of claim sizes (so that a claim of size X\ is 
made at time To, etc). We assume that the random variables from the above 
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3<N(t) 




N(t) — inf{j > : T + ■ ■ ■ + 2} > t}, 



sequences are jointly independent, with {Tj}j>\ and {Xj}j>i being i.i.d. se- 
quences. The goal of the present note is to derive an explicit formula for the 
distribution of the ruin time 

r = inf{t > : U(t) < 0} 

in the special case when the X^s follow the exponential distribution. 

When claims occur according to a Poisson process and the claim size 
distribution is exponential, a solution for the distribution of the ruin time 
r has been known for many years. See, for example, [2], [H] and [9] for 
different solutions to this problem. In the recent paper [7j, the authors used 
analytical techniques to obtain an explicit formula for the density of r in 
the case when the Tj's have an Erlang distribution. In the present note, we 
present an alternative probabilistic method, which enables one to derive such 
an explicit formula in the more general case when the Tj's follow an arbitrary 
distribution. 



2 The main result 

We assume that the claim sizes Xj follow the exponential distribution with 
parameter A > 0: 

F(Xj >x) = e~ Xx , x > 0, (1) 

while the positive random variables To and Xi(= Tj, j > 1) have densi- 
ties fo(t) and f(t), respectively. By g * h we denote the convolution of the 
functions g, h defined on (0, oo): 

( g * h)(t) = / g(t-v)h(v)dv, 
Jo 

and by g* n = g*( n_1 ) * g, n > 2, the n-fold convolution of g with itself. 

Theorem 1 Under the above assumptions, the ruin time r has a (defective) 
density p T (t) given by 



p T (t) = e~ x ^ { f (t) 

E A " ( ^, ctr " 1 Hr * / )(t) + <r * h)m }, (2) 



n=l 

where fi(t) = tf (t). 
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Proof. The idea of the proof is similar to the one used in [1] : first we will 
translate our problem into the problem of the crossing of a linear boundary 
by the pure jump process U°(t) = U(t) — ct and then swap the roles of the 
time and space coordinates. Then we notice that the generalised inverse of 
the function U°(t) is nothing else but the trajectory of a compound Poisson 
process. Eventually, the original problem proves to be equivalent to finding 
the distribution of the hitting time of a level by a skip-free Levy process, of 
which the solution is well-known and is given by Kendall's identity (see e.g. 
§ 12, Theorem 1 in [3], or [5]). 

(i) We will assume in parts (i)-(ii) of the proof that T = v = const (which 
is equivalent to conditioning on T , but is more convenient from a notation 
viewpoint). 

As we have just said, it is easily seen that, for the pure jump process 
U°(t) = U(t) -ct = u- J2 X i' 

j<N(t) 

one has 

t = inf{* > : U°(t) - (-ct) < 0}. 

Next we 'translate' the origin to the point (v, u) and swap the roles of co- 
ordinates by introducing the new 'time' s = u — x and 'space' y = t — v 
(where t and x respectively represent the original time and space). In the 
new system of coordinates, the trajectory of our process {U°(t)} is again a 
pure step function, which starts at zero at 'time' s = and has jumps of 
sizes Ti,T 2 , T 3 , . . . , at 'times' X\, X\ + X 2 , X\ + X 2 + X 3 , .... Due to our 
assumption (JT1), this will be a trajectory of the compound Poisson process 

Z°(s)= Yl T ^ 

k<M(s) 

where 

M(s) = inf{/c > 1 : X x + ■ ■ ■ + X k > s} - 1 

is a Poisson process with rate A. The distribution of the r.v. Z°(s) with s > 
has an atom e~ Xs at zero and a density on (0, oo) given by 

Pz°( S} (y) = e- Xs J2-^-r n (y)> y> Q - ( 3 ) 

71=1 
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To a crossing of the (lower) linear boundary x = —ct by the process 
{U°(t)} at time r (this necessarily is a jump epoch) there corresponds a 
(continuous) crossing of the (again lower linear) boundary 

y = s/c — (v + u/c), s > 0, 

by the process {Z°(s)} at 'time' a = u + cr, so that 

T = (a - u)/c. (4) 

Finally, we notice that a is the crossing time of the (lower) level —(v+u/c) 
by the process Z(s) = Z°(s) — s/c, which is clearly a skip-free in the negative 
direction Levy process. 

Figure 1 illustrates the translation of the original problem. The original 
surplus process starts at level u — 1, and ruin occurs at the fifth claim. Ro- 
tating the figure anti-clockwise through 90 degrees we see the corresponding 
path of the pure jump process {Z°(s)}. 




Figure 1: Original and translated processes 



(ii) Therefore, provided that Z(s) has a density Pz(s)(y) at the point 
y = —(v + u/c), the crossing 'time' a also has a density p a (s) at the point s, 
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which is given by Kendall's identity (see e.g. § 12, Theorem 1 in [3], or [5]): 

v + u/c 



: PZ(.s)(-(v + U/C)). 



This together with ^ implies that, for s > u + cv, the r.v. a has the density 



71=1 



??■! 



Therefore, it follows now from (HI) that, given To = v, for £ > v the stopping 
time t has a conditional density given by 



... . . u + cv 

Pr{t\v) = cp a {u + ct) = — — e 

u + ct 



Hu+ ct) y (X(u + ct)) n _ 



n=l 



(hi) To obtain the density of r in the general case, we observe that r > T 
always and so, using Q, 

P(r < t) = P(T = r < £) + P(T < r < t) 

?(u + cv-X 1 < 0)f (v)dv + / F(v < t < t\ T = v)f (v)dv 



- x(u+cv) fo(v)dv + 



p T (r\v)dr 



fo(v)dv. 



Differentiating both sides and substituting the representation for p T (t\v) from 
([5]) yields the density of r: 



Pr (t) = e^(" +rf Vo(t) + / p T (t\v)f (v)dv 



-\(u+ct) 



oo 



(A(« + ct)) 



71 /•* 



7t=l 



(u + cv)f* n (t-v)f (v)dv 



(the change of the order of integration/summation is justified as the integrand 
is a non- negative function). As the last expression is equivalent to the RHS 
of (J2|, the theorem is proved. □ 
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3 Examples 



3.1 Gamma inter-claim times 

Let us first consider the situation where claims occur according to an ordinary 
renewal process, so that each Tj, j = 0, 1,2, ... has density function 

orij.n-1 -f3t 

where n > and (3 > 0. It is well known that 

/On(m+lW(m+l)-l -fit 
f*(m+l)(f\ _ r l _ c 

J W r(n(m + l)) 

and it is straightforward to show that 

„ an(m+l)+lj.n(m+l) p - (3t 

J Jl{ ! p r(n(m+l) + l) ' 
Then formula ([2]) gives 

P (t) fp t \n-l U P e ~ KU+Ct) ~ Pt ^ \ m (u + Ct) m (pt) nm 



u + ct ^— ' m! r(n(m + 1)) 

m=0 



cne- A ( M+c *)- /3 * ^ A m (w + ct) m (/ft) 



u + ct ^— ' m! r(n(rn + 1) + 1") 

m=0 



In the special case when n is a positive integer, we can compute this as 

^ e -A( M+ ct)-ffl (/?t)n -i / / 1 n - 1 X(u + ct){(3ty 

Pr{t) = , 7 rr- Mofn M + -, .",H 



u + ct V V ' n n n n 

+ ct ^ 1 + -1 + -,-, 1 + -;- , 6 



where 



,F q ( Bl , B 2 , ...,B p , Cl ,C 2 ,... C q , Z) - ^ {Cl)m{C2)mACq)m ^ 
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is the generalised hypergeometric function (and (a) n = T(a + n)/T(a) is 
Pochhammer's symbol). Formula (Is]) follows from the identity 



r((n(m + l) + l) n nm JJ[r(m + l + ^)' 

which can be derived by applying the multiplication formula of Gauss as 
described in [7j. 

Formula (|6]) is in a different form to the formula for p T (t) derived in [7] . 
A comparison of these two formulae for p T (t) yields the identity 



1 n-1 \(u + ct)(pty 

O-Tn J-, J- H , -L + 



n n n' 



12 n A(« + ct)(/%)' 

— o-^n J- H j J- H > A + 



n n n n' 



Ku + ctmr J?L oF J 2 + l, 2 + l 2 + "A(« + W 



(2n)! \ n' n' n' n r 



In the special case n — 1, by writing z = a/4A/3£(u + ct) this identity reduces 
to the well-known result (e.g. pQ) 

h{z) - -h(z) = h{z), 
z 

where I v is the modified Bessel function of order v. 

Next, let us consider the special case when n = 2, and let us further 
assume that claims occur according to a stationary renewal process, so that 
the distribution of To is the equilibrium distribution of T\. Then we find that 

fo(t) = \e-^ (1 + (3t) = Itfe-* + pHe'V), 

giving 

-y / p2m+1^2m e -l3t p2m+2^2m+l e -f3t 

f m * fo(t) = ^ [ r(2m + l) + r(2m + 2) 

Further, 

f 1 (t) = tf (t) = ^te^ t + p 2 t 2 e^ t ), 
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giving 

^ p2m+lj-2m+l e -f3t j^2m+2^2m+2 e ~f3t 

f m * fi{t) = 2 r(2m + 2) + T(2m + 3) ' 
Then formula ^ gives 



Pt{ > " l /0 ^ + u + ct Z_, m! 2 \ T(2m+1) + T(2m + 2) 



m=l 

oo 



^o. + r+l-^ m) \ 9 rT9m. -I- 9"l + 



u + ct^ ml \2 T(2m + 2) T(2m + 3) 

m=l 



and we can incorporate fo(t) into the sums so that both start at m — 
0. For computational purposes we can write this in terms of generalised 
hypergeometric functions as 

Pr(*) = o/ , « 0^2 ^, 1; 



2(u + ct) V V 2 ' ' 4 



(7) 



Table [T] shows some values of finite time ruin probabilities when A = 1, 
(3 — 2 and c = 1.1. We use the notation ip(u,t) to denote the probability 
of ruin by time t from initial surplus u when the density of r is given by 
formula (|6| with n — 2, and ip e (u,t) denotes the corresponding probability 
when the density of r is given by formula ([7]). These values have been found 
by integrating the density functions numerically using Mathematica. We 
can observe from this table that for each combination of u and t, the finite 
time ruin probability is greater when the distribution of T is the equilibrium 
distribution of 7\. This arises because both the mean and variance of T are 
smaller than the corresponding values for T\. 



3.2 Mixed exponential inter-claim times 

Let us now consider the situation when the distribution of each Tj, j = 
0, 1, 2, is mixed exponential with density function 

f(t) = f (t)=pae- at + q{3e-P\ 
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t 




MM 


^(10, t) 


M^-M 


V(20,t) 


M 2 M 


20 


0.7973 


0.8463 


0.0457 


0.0509 


0.0009 


0.0010 


40 


0.8332 


0.8735 


0.1008 


0.1082 


0.0060 


0.0066 


60 


0.8481 


0.8848 


0.1387 


0.1469 


0.0138 


0.0148 


80 


0.8564 


0.8912 


0.1651 


0.1737 


0.0218 


0.0232 


100 


0.8618 


0.8952 


0.1842 


0.1930 


0.0292 


0.0309 



Table 1: Finite time ruin probabilities. 



where < p < 1, q — 1 — p, and > a > 0. Following ideas in [10], it is 
shown in [6] that the m-fold convolution of / with itself as can be written as 



/*™(t)=^ 7mje ( m + j,/3;t), 
3=0 

where e(m,/3;t) denotes the Erlang(m) density with scale parameter j3 and 



We can find a similar type of expression for f* m * fi(t) by using Laplace 
transforms. For a function w, let 



Then 



and 



w(s) 



e~ st w(t)dt. 



J[S) a+s + 0+s 



fi(s) 



pa 



+ 



q0 



leading to 

" _ ~\m _ 

f(s) fi(s) 



P 

a 



r=0 



p r q m-r 



+-Y 



r=0 



p r q m-r 



a + s 
a 



r+2 



a + s 







13 + s 







+ s 



m—r+2 



Hence 



Jo T(r + 2) r(m-r) 



a JL / m \ r t a r v r ~ 1 e~ ay 3 m ~ r+2 (t — y) m - r+1 e ~ f3( - t ~ y ^ 

+^E( WW y v( \ - — I iw — d f 

f3f^ Q \ r J Jo r ( r ) r(m-r + 2) 

7^V r / r(m + 2) 

f^o \ r J r(m + 2) 

If we now replace the i-F\ functions by their series representations, we find 
after a small amount of manipulation, that 

oo 

/*- * f t (t) = J2 Vi,me(™ + * + 2, /3; t), 

i=0 

where 

^ = ^ (1 " E (7) ^ (g) r ((r + 2)(r + l)ap + /3g) . 



Thus, we have formulae for all the ingredients in formula (pq). 

Figure 2 shows three plots of the density of r when u = 10, c = 1.1 and 
the parameters of the mixed exponential distribution are as in Table [2] In 
Figure 2, Case A is illustrated by the dotted line, Case B by the solid line, 
and Case C by the bold line. We observe that the ordering of these three 
plots leads to highest finite time ruin probabilities for Case A and lowest for 
Case C, consistent with the ordering of the three values of V [To]. 

Case a f3 p E [T ] V[T ] 



A 2/5 2 1/4 1 5/2 
B 1/2 2 1/3 1 2 
C 3/5 2 3/7 1 5/3 

Table 2: Parameters of mixed exponential distributions. 
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o.oos 



Q.Q04 



Figure 2: Densities of time to ruin for mixed exponential inter-claim times. 
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